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ABSTRACT 

A population of the aquatic gastropod Bythinella that does not conform to any 
described hypogean or epigean species of this genus was discovered in a cave 
stream in southern France (Folatiere cave). Members of this population lacked 
body pigments and eyespots, which suggests adaptation to a subterranean 
environment. Mitochondrial COI sequences revealed two lineages within this 
population, separated by a genetic divergence higher than the inter-specific 
threshold commonly admitted for hydrobioid species. In contrast, a phylogeny 
based on the nuclear ITS 1 gene clustered the sequenced hypogean specimens 
in a single clade. Moreover multivariate analyses performed on shell parameters 
did not enable discrimination between hypogean specimens attributed to each 
of these two mitochondrial lineages. The high COI variability within the same 
population and conflict between mtDNA and nDNA trees could be explained 
by a scenario of introgression by hybridization. Consequently, based on nuclear 
gene phylogeny, shell and anatomical characters, and ecological evidences, we 
hypothesize that the easily recognizable hypogean morph of the Folatiere cave 
can be attributed to a distinct and new species - here described and named 
B. navacellensis n. sp. Our extensive sampling suggests that B. navacellensis 
n. sp. is distributed within a small delimited and cohesive area located around 
southern Larzac, covering c. 300 km 2 . Considering its distribution, restricted 
to a few subterranean aquifers only, B. navacellensis n. sp. can be categorised as 
a vulnerable species, VU D2. 
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RESUME 

Relations phylogenetiques et description d’une nouvelle espece stygobie du genre 
Bythinella (Mollusca, Gastropoda, Caenogastropoda, Amnicolidae) du sud de la 
France. 

Une population de gasteropodes aquatiques du genre Bythinella , ne correspondant 
a aucune espece decrite, a ete decouverte dans une riviere souterraine dans le 
sud de la France (la grotte de la Folatiere). Les individus de cette population 
sont depigmentes et anophtalmes, suggerant une adaptation a l’environnement 
souterrain. Les sequences du gene mitochondrial COI ont revele deux lignees au 
sein de cette population, avec une divergence genetique superieure a la moyenne 
de variability interspecifique communement admise pour les hydrobies. Toutefois, 
une phylogenie sur la base du gene nucleaire ITS 1 regroupe tous les individus 
sequences dans le meme clade. De plus, les analyses multivariees realisees sur 
les parametres morphologiques de la coquille n’ont pas permis de discriminer 
ces deux lignees mitochondriales. La forte variability du gene COI au sein de 
la meme population et le conflit entre les arbres de l’ADN mitochondrial et 
nucleaire pourraient s’expliquer par un scenario d’introgression par hybridation. 
En consequence, en nous appuyant sur la phylogenie du gene nucleaire, les 
caracteres anatomiques et morphologiques de la coquille et les particularites 
ecologiques, nous posons fhypothese que ce morphe facilement reconnaissable 
peut etre attribue a une espece nouvelle - decrite ci-apres et nominee Bythinella 
navacellensis n. sp. Notre echantillonnage intensif suggere que B. navacellensis 
n. sp. est limitee a une petite zone cohesive des causses du Larzac meridional et 
de Blandas, couvrant environ 300 km 2 . En raison de sa distribution restreinte 
a quelques aquiferes souterrains, B. navacellensis n. sp. peut etre categorisee 
comme vulnerable YU D2. 


INTRODUCTION 

The genus Bythinella Moquin-Tandon, 1856 consists 
of small hydrobioid springsnails widespread from 
northern Spain to eastern Turkey. This group is 
considered as one of the most diverse amongst 
European hydrobioids, with 91 species and subspecies 
(Boeters 2006, 2008; Gloer 2006; Gloer & Pesic 
2006; Fiaase et al. 2007; Boeters & Falkner 2008; 
Fauna Europaea Web Service 2009), and France, 
with its 49 alleged valid taxa, is its main centre of 
species richness. 

Bythinella species live mainly in springs and in 
the uppermost courses of small streams (Falniowski 
1987). Fiowever, five taxa from France are stygobite, 
i.e. living strictly in underground waters: B. pupoides 
phreaticola Bernasconi, 1989 and B. geisserti Boeters & 


Falkner, 2003 were described from hyporheic zones 
in the French Jura region; B. padiraci Locard, 1903 
was described from the Padirac karstic network in 
the Lot (see Bichain etal. 2004); B. bouloti Girardi, 
Bichain & Wienin, 2002 and B. galerae Girardi, 
Bichain & Wienin, 2002 are both restricted to a 
few caves in the Grands Causses of Lozere. 

Because access to, and sampling in, hypogean 
habitats are difficult, Bythinella species have been 
described mainly from very few living animals or 
from empty shells only sometimes collected outside 
the subterranean networks. Consequently species 
delimitations and/or ecological status (stygobite 
vs. stygoxene, i.e. strict vs. occasional subterranean 
species) of hypogean Bythinella populations are 
still unclear and subject to various opinions. For 
example, Bertrand (2004) recognized five stygobite 
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Bythinella taxa whereas Falkner {in Ferreira et al. 
2003) recognized seven. Additionally, Bernasconi 
(2000) recognized Bythinella eutrepha (Paladilhe, 
1867) as an occasional hypogean species, widespread 
in southeastern France whereas Falkner {in Ferreira 
etal. 2003) considers it as a strictly stygobite species 
endemic to the Lez catchment. 

One of us (VP, August 2003) collected c. 45 living 
Bythinella specimens in the underground waters of 
the Folatiere cave (Fig. 1A, B) in southern France. 
Morphologically, they did not correspond to any 
described species. Subsequently, empty shells of 
this morph were also sampled in 21 other subter¬ 
ranean localities of the Larzac and Blandas plateaus, 
a large karstic zone located north of Montpellier 
(France). No living specimens were found in any 
of these localities. 

In order to explore the taxonomic status of this 
morph, five specimens collected in the Folatiere cave 
were first included in a molecular phylogenetic study 
that focused on species delimitation problems within 
the genus Bythinella (Bichain etal 2007a). The results 
highlighted two divergent mitochondrial lineages 
among syntopic specimens sampled in the cave. 
Four specimens clustered with specimens unambigu¬ 
ously attributed to the epigean species B. eurystoma 
(Paladilhe, 1870) and one specimen constituted a 
distinct and unnamed lineage. No shell characters 
allowed morphological discrimination among these 
syntopic hypogean specimens that would support 
this mitochondrial lineage divergence. 

In contrast, the nuclear ITS1 sequences obtained 
from two hypogean specimens were deeply diver¬ 
gent from those of epigean specimens attributed to 
B. eurystoma. ITS1 thus indicated that B. eurystoma 
would be a species restricted to epigean habitats 
whereas all hypogean specimens from Folatiere cave 
would belong to a distinct stygobite species. 

The few individuals analyzed so far from the 
Folatiere cave, coupled with the fact that substitution 
rate is higher in the mitochondrial genome than in 
the nuclear one, may suggest that the genetic break 
observed by Bichain etal. (2007a) could be due to a 
sampling effect: two divergent haplotypes may have 
been randomly sampled in the population. 

Here, we increased the dataset of Bichain et 
al. (2007a) by including sequences from 16 ad¬ 


ditional specimens sampled in the Folatiere cave. 
Our objectives were 1) to test if the two divergent 
mitochondrial lineages still persist when further 
individuals are analyzed, and 2) to examine diver¬ 
gences/similarities in shell morphology between the 
epigean and hypogean populations and within the 
hypogean population. 

Abbreviations 

a.s.l. above sea level; 

BI Bayesian Inference; 

COI cytochrome oxydase subunit I; 

DFA Discriminant Function Analysis; 

ITS1 Internal Transcriber Spacer 1; 

mtDNA mitochondrial DNA; 
nDNA nuclear DNA; 

PCA Principal Component Analysis; 

PC Principal Component; 

L lens; 

L a aperture height; 

L s h shell height; 

L} w last whorl height; 

R retina; 

W a aperture width; 

W LW last whorl width; 

W s h shell width; 

WN whorls number; 

WL Wilks’ lambda; 

MNHN Museum national d’Histoire naturelle, Paris. 


MATERIAL AND METHODS 

Study area and sampling methods 
The study area is located in the Lez and Herault 
catchments, southern France (Fig. 1). In order to 
collect hypogean Bythinella , sampling was performed 
in caves, springs and in the hyporheic zone. For 
each of these biotopes, we used several collecting 
methods: 1) in caves, living animals were collected on 
the spot by sight and sediment was sampled; in the 
Folatiere cave, divers collected sediment down to a 
depth of 22 m; 2) in springs, sediment was collected 
and drifting fauna was sampled using drifting nets 
during flood events; and 3) the hyporheic zone was 
investigated with the Bou-Rouch pumping method 
(Bou & Rouch 1967). Cave and spring sediments 
were first soaked in tap water and left in the dark 
for a few days in order to leave enough time for 
live animals to reach the surface of the sediment. 
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Fig. 1 . — Study area and sampling sites: A, view of the Folatiere cave entrance; B, view inside the Folatiere cave (type locality of 
Bythinella navacellensis n. sp.) showing the typical erosion resulting from the active network; C, location of sampling sites (station 
code and details on localities are given in Table 1); Thick lines, main rivers; o, sampling sites (sites with occurrence of B. navacellensis 
n. sp. are numbered, data from rivers alluvium are not shown); ▲, sampling sites of B. eurystoma (Paladilhe, 1870). The geographical 
range of B. navacellensis n. sp., inferred from hydro-geological data (see text), is darkened. 
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Table 1. — List of Bythinella navacellensis n. sp. and B. eurystoma (Paladilhe, 1870) sampling sites. Site numbers from 1 to 21 and 
26 to 28 are plotted on Figure 1. 


Sites Species Sample sites description Latitude Longitude Hydrogeological zone 

no. 


1 

B. 

navacellensis 

n. 

sp. 

Springs of Lasfons (1) 

696 400 

1 

885 

700 

Causse de 

Blandas (north) 

2 

B. 

navacellensis 

n. 

sp. 

Springs of Lasfons (2) 

697 800 

1 

885 

900 

Causse de 

Blandas (north) 

3 

B. 

navacellensis 

n. 

sp. 

Folatiere cave (locality 382) 

695 750 

1 

874 

280 

Causse de 

Blandas (south) 

4 

B. 

navacellensis 

n. 

sp. 

Spring of Folatiere cave 

695 652 

1 

874 

182 

Causse de 

Blandas (south) 

5 

B. 

navacellensis 

n. 

sp. 

Gourneyras cave 

695 600 

1 

874 

110 

Causse du 
(north) 

Larzac meridional 

6 

B. 

navacellensis 

n. 

sp. 

Gourneyrou cave 

695 720 

1 

874 

040 

Causse du 
(north) 

Larzac meridional 

7 

B. 

navacellensis 

n. 

sp. 

Small spring on the side Vis 

696 440 

1 

873 

240 

Causse du 

Larzac meridional 






river 





(north) 


8 

B. 

navacellensis 

n. 

sp. 

Magnanerie cave 

697 830 

1 

873 

010 

Causse de 

Blandas (south) 

9 

B. 

navacellensis 

n. 

sp. 

La Tuilede cave 

698 180 

1 

873 

510 

Causse de 

Blandas (south) 

10 

B. 

navacellensis 

n. 

sp. 

Spring of Grenouillet 

697 550 

1 

872 

670 

Causse du 
(north) 

Larzac meridional 

11 

B. 

navacellensis 

n. 

sp. 

Springs of Mas-Neuf 

701 200 

1 

874 

800 

Causse de 

Blandas (south) 

12 

B. 

navacellensis 

n. 

sp. 

Calavon cave 

701 400 

1 

875 

400 

Causse de 

Blandas (south) 

13 

B. 

navacellensis 

n. 

sp. 

Spring in the Cirque du bout 
du monde 

684 700 

1 

865 

500 

Causse du 

Larzac (south) 

14 

B. 

navacellensis 

n. 

sp. 

Spring of Saint-Laurent 

684 000 

1 

865 

000 

Causse du 

Larzac (south) 

15 

B. 

navacellensis 

n. 

sp. 

Banquier cave 

685 100 

1 

864 

400 

Causse du 

Larzac (south) 

16 

B. 

navacellensis 

n. 

sp. 

Spring of Parlatge (1) 

687 100 

1 

864 

600 

Causse du 

Larzac (south) 

17 

B. 

navacellensis 

n. 

sp. 

Spring of Parlatge (2) 

687 300 

1 

864 

800 

Causse du 

Larzac (south) 

18 

B. 

navacellensis 

n. 

sp. 

Springs of Ruisseau de la 
Tour 

689 600 

1 

860 

600 

Causse du 

Larzac (south) 

19 

B. 

navacellensis 

n. 

sp. 

Springs of Valrousse 

688 800 

1 

861 

900 

Causse du 

Larzac (south) 

20 

B. 

navacellensis 

n. 

sp. 

Springs of Saint-Saturnin 
de Lucian 

690 200 

1 

859 

600 

Causse du 

Larzac (south) 

21 

B. 

navacellensis 

n. 

sp. 

Spring of Arboras 

692 800 

1 

857 

000 

Causse du 

Larzac (south-east) 

22 

B. 

navacellensis 

n. 

sp. 

Alluvium of the Vis river 

695 000 

1 

876 

900 

Unknown 


23 

B. 

navacellensis 

n. 

sp. 

Alluvium of the Vis river 

699 490 

1 

873 

230 

Unknown 


24 

B. 

navacellensis 

n. 

sp. 

Alluvium of the Vis river 

694 700 

1 

877 

800 

Unknown 


25 

B. 

navacellensis 

n. 

sp. 

Alluvium of the Lergue river 

679 880 

1 

862 

120 

Unknown 


26 

B. 

eurystoma 



Spring near St-Julien- 
de-la-Nef (locality 332) 

— 


— 


Unknown 


27 

B. 

eurystoma 



Spring near Montdardier 
(locality 333) 

— 


— 


Causse de 

Blandas (north) 

28 

B. 

eurystoma 



Spring near Saint-Maurice- 

- 


- 


Causse de 

Blandas (south) 


Navacelles (locality 329) 


After collecting the live specimens, the sediment was 
dried and poured into water again, thus allowing 
empty shells to surface. Sampling involved more 
than 60 localities of which only one revealed living 
hypogean animals (Folatiere cave, Fig. 1A, B) and 
21 (from caves and springs only) revealed empty 
shells attributable to the B. navacellensis n. sp. 
morphotype (Fig. 1C; Table 1). All live specimens 
were fixed in 95% ethanol. 

Sixteen living animals sampled from the Folatiere 
cave were analyzed both genetically and morpho¬ 


metrically. For the latter analyses, we included 
empty shells from the Folatiere cave and specimens 
of B. eurystoma from three epigean localities in the 
Larzac region (Fig. 1C; Table 1). 

Molecular analyses 

Total DNA was extracted from whole individuals (shell 
included) using the QLAGEN DNeasy kit (Qiagen Inc., 
Hilden, Germany). Partial COI mtDNA was amplified 
using the universal primers H2198 and LI490 (Folmer 
etal. 1994). PCR reactions were performed in a final 
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volume of 25 pi, using c. 2.5 ng of template DNA, 
1.25 mM MgCl 2 , 0.3 pM of each primer, 0.13 mM 
of each nucleotide, 5% DMSO and 1.5 unit of Taq 
polymerase (MP Biomedical Qbiogene Inc., Illkirch, 
France). Amplification products were generated by an 
initial denaturation step of 4 min at 94°C, followed 
by 35 cycles at 94°C for 30 s (denaturation), 52°C for 
40 s (annealing) and 72°C for 40 s (extension), and 
a final extension of 10 min at 72°C. PCR products 
were purified using the QIAquick PCR purification 
kit (Qiagen Inc.). Sequencing was performed by 
Genoscreen Inc. (Lille, France) on an ABI prism 
3730XL (Applied Biosystems, Foster City, California) 
in both directions to confirm authenticity of each 
sequence. 

We added 16 COI sequences (GenBank num¬ 
bers EF372562-EF372577) to the global dataset of 
Bichain etal (2007a) (92 sequences, GenBank num¬ 
bers EF016165-EF016248). Three amnicolid COI 
sequences from GenBank (AF367652, AF322409, 
AF213348) were used to root our analyses. COI 
sequences were aligned unambiguously by eye using 
BioEdit 5.0.0 (Hall 1999). No stop codon or shifts 
in the reading frame were detected when translated 
into amino acids. 

Phylogenetic analyses were conducted using Baye¬ 
sian inference (BI) as implemented in MrBayes v3.1.1 
(Ronquist & Huelsenbeck 2003). Using the AIC 
criterion, the GTR + I + T model (Yang 1994) was 
selected by MrModeltest 2.2 (Nylander 2004); we 
subsequently used this model for all the analyses. 
Flat priors were used for all parameters. Four Markov 
chain Monte Carlo simulations (one heated, three 
cold) were run for 2 million generations and pa¬ 
rameters were sampled every 100 generations. The 
first 2 x 10 5 generations (2000 trees) were discarded 
(“burn-in” period) and the posterior probabilities were 
estimated from the remaining sampled generations. 
We checked that the potential scale reduction factor 
approached 1.0 for all parameters. We used Tracer 
vl .4 (Rambaut & Drummond 2007) to ascertain 
that our sampling of the posterior distribution has 
reached enough effective sample size. 

Morphological analyses 

Shells were placed on an adhesive support in a 
standard position (Bichain et al 2007b) and then 


digitalized with a graduated scale using a stereo¬ 
microscope connected to a digital camera. Seven 
parameters were recorded on each picture using 
ImageTool 3.00 (Wilcox et al 2002): height and 
width of the shell, height and width of the last whorl, 
height and width of the aperture, the whorls number 
visible on each picture. Four ratios (W^/L^, L[ w / 
L s h, L a /L s h, W a /L a ) commonly used in the alpha- 
taxonomy of hydrobioid spring-snails (Bernasconi 
2000; Haase et al. 2005) were calculated. 

We performed multivariate analyses on these 11 
shell parameters measured on: 1) 72 epigean shells 
attributed to B. eurystoma ; 2) the 16 new specimens 
from the Folatiere cave involved in our molecular 
study, which were digitized before the DNA extrac¬ 
tions and; 3) 19 empty shells from this cave that 
could not be specifically attributed. 

We first performed a Principal Component Analy¬ 
sis (PCA) to explore the distribution of the shell pa¬ 
rameters between epigean and hypogean specimens. 
Then Discriminant Function Analyses (DFA) were 
used to test statistically: 1) the differences between 
hypogean shells attributed genetically to divergent 
haplotypes, and 2) the differences between hypo¬ 
gean and epigean individuals. Wilks’ lambda used 
in an ANOVA (F) test of mean differences was 
used to test if the discriminant model as a whole 
was significant. In such cases, each variable was 
then tested using WL to determine which vari¬ 
able differed significantly in mean value between 
discriminated groups. 

All analyses were performed with Statistica® 6.0 
and the level of significance used was set to 5%. 

RESULTS 

Phylogenetic analyses 

The Bayesian tree obtained from 111 COI sequences 
agreed with the phylogenetic results of Bichain etal 
(2007a) and confirmed the strong genetic divergences 
between the two hypogean mitochondrial lineages 
from the Folatiere cave. Eight new COI sequences 
clustered (Fig. 2A) in the clade A (. B . navacellensis 
n. sp. and B. eurystoma) and eight others in the 
clade B (. B . navacellensis n. sp. and B. cebennensis 
(Dupuy, 1849)). 
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j- Aniane (339-4) 
H—Cabrier (338-2) 
|j- Brissac (334-4) 
p^-Brissac (334-2) 
rAniane (339-1) 
Va Bueges (330-2) 
-La Bueges (330-5) 
-La Bueges (330-1) 

- Brissac (334-5) 

L St-Laurent (331-2) 

- Cabrier (338-3) 

-St-Laurent (331-1) 



B. cebennensis 


Clade B 


r- Folatiere (382-D) 
P- Folatiere (382-4) 
i—[Folatiere (382-A) 

J * Folatiere (382-E) 

|j-Folatiere (382C) 

I - Folatiere (382-F) 

- Folatiere (382-G) 
Folatiere (382-H) 

L Folatiere (382-B) 


Clade A 


j—Folatiere (382-2) 

I"*- Folatiere (382-1) 
n_ Folatiere (382-J) 

L Folatiere (382-N) 

Lr Folatiere (382-L) 

^ Folatiere (382-M) 

- Folatiere (3821) 
p Folatiere (382-1) 

I— Folatiere (382-P) 
T-Folatiere (3823) 
^-Folatiere (382-0) 

— Folatiere (382-K) 

I—Montdardier (333-1) 

J-St-Julien-Nef (332-4) 

L Montdardier (333-3) 

I—St-Julien-Nef (332-5) 
St-Julien-Nef (332-3) 
Navacelles (329-5) 
Navacelles (329-1) 

Navacelles (329-2) 

_p Montdardier (333-5) 

~Montdardier (333-4) 


_T st 

jLn 

l-Na\ 
I— Ma 



B. eurystoma 


Fig. 2. — Bayesian phylogenetic tree of Bythinella specimens based on 111 COI sequences: A, detail of the COI global Bayesian 
tree implying the clade A (B. navacellensis n. sp. and B. eurystoma (Paladilhe, 1870)) and clade B (B. cebennensis (Dupuy, 1849) and 
B. navacellensis n. sp.) (shell outlines represented). Name and code of each locality are indicated in Table 1 for B. eurystoma and 
B. navacellensis n. sp. Thick bars, exclusively hypogean species. Values of branch robustness are indicated for each node; model 
parameters for Bl analyses: I = 0.228; T 4 = 0.124; [AC] = 0.035; [AG] = 0.459; [AT] = 0.050; [CG] = 0.048; [CT] = 0.291; [GT| = 0.026; A = 
0.327; C = 0.221; G = 0.115; T = 0.327; B, detail of the ITS1 global Bayesian tree modified from Bichain etal. (2007a). The individuals 
of B. navacellensis n. sp. for which we obtained both COI and ITS1 sequences are indicated in bold. Values of branch robustness are 
indicated for each node. 


Within the clade A, there were two haplotypes 
diverging by four substitutions. Each haplotype corre¬ 
sponded respectively either to the epigean (. B . eurystoma 
lineage, 22 sequences) or hypogean (. B . navacellensis 
n. sp. lineage, 12 sequences) individuals. Within the 
B. navacellensis n. sp. lineage belonging to the clade 
B (nine sequences), three haplotypes diverged by 
two substitutions at most. The intra-clade p-distance 
ranged from 0.1% (within clade A) to 0.4% (within 
B. navacellensis n. sp. lineage from clade B) and the 
inter-clade p-distance [clade A vs. B. navacellensis 
n. sp. lineage from clade B] reached 4.7%. 


Multivariate analyses 

Our global dataset included 11 parameters measured 
on 107 shells. From epigean localities, 72 shells 
were attributed to B. eurystoma, , a species genetically 
and morphologically distinct from B. cebennensis. 
From the Folatiere cave 16 were attributed to 
B. navacellensis n. sp. and 19 empty shells were 
not specifically attributed. Basic statistics (mean 
and standard error) of the 11 shell parameters are 
given in the Table 2. 

To explore the pattern of shell variation, a PCA 
was carried out on the entire dataset (107 epigean 
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Table 2. — Basic statistics of 11 shell parameters fo rBythinella eurystoma (Paladilhe, 1870) and B. navacellensis n. sp. Measurements 
in mm. 


Bythinella eurystoma 


N = 72 

L sh 

L| W 

La 

W lw 

W sh 

W a 

W sh /L sh 

l-a/l-sh 

l-lw/l-sh 

W a /L a 

WN 

mean 

2.666 

1.906 

1.210 

1.292 

1.574 

0.999 

0.592 

0.455 

0.717 

0.827 

3.361 

Minimum 

2.269 

1.678 

1.064 

1.106 

1.379 

0.904 

0.537 

0.409 

0.636 

0.726 

3.000 

Maximum 

3.108 

2.141 

1.366 

1.476 

1.812 

1.119 

0.665 

0.489 

0.780 

0.918 

4.000 

Standard error 

0.198 

0.102 

0.075 

0.073 

0.091 

0.050 

0.034 

0.018 

0.029 

0.039 

0.319 

Specimens from the Folatiere cave [16 living animals and 19 empty shells attributed to B. navacellensis n. sp.] 

N = 35 

Lsh 

L| W 

La 

W lw 

W sh 

W a 

W sh /L sh 

La/Lsh 

Liw/L S h 

W a /L a 

WN 

mean 

2.460 

1.703 

1.038 

1.350 

1.591 

0.867 

0.650 

0.423 

0.695 

0.835 

3.689 

Minimum 

2.038 

1.489 

0.864 

1.193 

1.335 

0.675 

0.575 

0.389 

0.646 

0.762 

3.500 

Maximum 

2.960 

1.961 

1.209 

1.591 

1.914 

1.057 

0.758 

0.475 

0.740 

0.909 

4.000 

Standard error 

0.272 

0.141 

0.092 

0.099 

0.126 

0.085 

0.041 

0.022 

0.026 

0.035 

0.246 


Specimens from the Folatiere cave [16 living animals attributed to B. navacellensis n. sp.] 


N = 16 

Lsh 

L| W 

La 

W lw 

W sh 

W a 

W sh /L sh 

La/Lsh 

L|w/L s h 

W a /L a 

WN 

mean 

2.613 

1.781 

1.085 

1.417 

1.648 

0.908 

0.633 

0.416 

0.684 

0.837 

3.844 

Minimum 

2.078 

1.520 

0.864 

1.275 

1.384 

0.719 

0.575 

0.389 

0.646 

0.766 

3.500 

Maximum 

2.960 

1.961 

1.209 

1.591 

1.914 

1.057 

0.691 

0.451 

0.732 

0.909 

4.000 

Standard error 

0.271 

0.141 

0.100 

0.093 

0.138 

0.096 

0.031 

0.016 

0.023 

0.041 

0.239 


and hypogean shells). The first (PCI) and second 
(PC2) principal components respectively account for 
44.9% and 27.4% of the total variation (Fig. 3A). 
Two groups of strongly correlated parameters were 
detected: 1) global size parameters (L^, W s ^, W^ 
Liw5 L a , W a ); and 2) ratios (L lw /L sh , L a /L sh ) repre¬ 
senting the relative size of the aperture and of the 
last whorl. The first set of variables was primarily 
associated with PCI. This first axis was thus inter¬ 
preted as a size factor. The second set of parameters 
was primarily associated with PC2. Additionally, 
W s h/L sh and W a /L a were weakly negatively corre¬ 
lated with PCI and the whorl number (WN) was 
strongly correlated with PC2. The projection of 
the individuals in the PCI x PC2 factorial plane 
revealed that B. eurystoma and all hypogean indi¬ 
viduals from Folatiere cave differed along PCI and 
PC2 (Fig. 3B). PCA failed to reveal any obvious 
difference between individuals of B. navacellensis 
n. sp. belonging to the two mitochondrial lineages 
revealed by the phylogenetic analyses. The PCA 
suggested that epigean B. eurystoma differed from 
all the hypogean specimens sampled in the Folatiere 
cave in having larger shells, larger last whorl, aper¬ 
ture different in shape and size, spire formed by a 
smaller number of whorls. 


The DFA carried out on the 16 living individuals 
from the Folatiere cave (eight attributed to clade A 
and eight attributed to clade B) and computed as 
“species” as the discriminant factor was not signifi¬ 
cant (WL = 0.144 and p = 0.239). The DFA car¬ 
ried out on the overall dataset and computed with 
“hypogean/epigean origin” (35 hypogean and 72 
epigean) as the discriminant factor was significant 
(WL = 0.063 andp < 0.05). Scores of WL with the 
probability associated with each variable indicated 
that both shell samples differed significantly in six 
parameters (L sh , L a , W sh , L a /L sh , L lw /L sh , WM). The 
overall reclassification score of the discriminating 
linear functions was 100%. 

The DFA carried out on the overall dataset, but 
excluding the 19 subterranean specimens not taxo- 
nomically attributed and computed with “species 
attribution” (72 B. eurystoma and 1 6 B. navacel¬ 
lensis n. sp.) as the discriminant factor was highly 
significant (WL = 0.569, p < 0.05). The overall 
reclassification score of the discriminating linear 
functions was 100%. 

The DFA results corroborate the PCA analysis. 
DFA indicated that morphological differences 
between hypogean and epigean shells were statisti¬ 
cally significant, whereas shell differences between 
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A B 



Fig. 3. — Multivariate analyses on shell parameters of Bythinella eurystoma (Paladilhe, 1870) and B. navacellensis n. sp.: A, correlation 
of the 11 parameters measured on the shell to the first two principal components, PCI and PC2; B, projection of the individuals in the 
first factorial plane PCI x PC2 (O, individuals of B. navacellensis n. sp. [from clade A]; □, individuals of B. navacellensis n. sp. [from 
clade B]; •, individuals of B. eurystoma: ■, undetermined individuals sampled in the Folatiere cave. All the hypogean individuals were 
sampled in the Folatiere cave, type locality of B. navacellensis n. sp.). 


hypogean mitochondrial lineages of B. navacellensis 
n. sp. were not. 

Hypogean individuals of the Folatiere cave (Fig. 4) 
seem to differ from epigean ones in having smaller 
shells and a more stocky shell shape, expressed 
by the lower values of L^/L^ and L a /L s h and the 
greater value of 

Species delimitation 

Our approach to species delimitation was based 
on a combination of morphological characters, 
ecological traits and genetic data. The population 
of the Folatiere cave exhibits a markedly peculiar 
shell and a clear ecological specificity. Additionally, 
previously published phylogenetic reconstruction 
based on the ITS 1 nuclear gene (Fig. 2B) suggests that 
this population constitutes a monophyletic lineage 
(Bichain etal. 2007a). The paraphyly suggested by 
the COI gene is interpreted as a consequence of 
introgression by hybridization (see Discussion). 
Based on nDNA tree, shell characters, anatomical 
features and ecological traits, we therefore consider 
this population as belonging to a distinct and new 
species named B. navacellensis n. sp. 


SYSTEMATICS 

Family AMNICOLIDAE Tryon, 1863 
Genus Bythinella Moquin-Tandon, 1856 

Bythinella navacellensis n. sp. 

Bythinella sp. - Bichain etal. 2007a: figs 2, 3, clade F + G 
(and Bythinella sp. 2 on the figure 5). 

Bythinella sp. - Prie 2008: fig. 1H. 

HOLOTYPE. — France, Folatiere cave, VII. 2003, V. Prie 
collector, 1 specimen (MNHN-9315; Fig. 5A). 

PARATYPES. — Same data as holotype, 19 specimens 
(MNHN-9316). 

Type locality. — France, Grotte de la Folatiere (Folatiere 
cave), town of Rogues, Gard departement, 3.5264E, 
43.8636N. 

ETYMOLOGY. — The specific epithet refers to the grandiose 
Cirque de Navacelles located near the Folatiere cave. 

VERNACULAR NAMES. — Bythinelle de Navacelles (French), 
Navacelles’ cave snail (English). 

Other material examined. — See Table 1. 
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Fig. 4. — Infra- and inter-specific shell variability of Bythinella eurystoma (Paladilhe, 1870) and B. navacellensis n. sp.: A-H, individuals 
sampled in the Folatiere cave and attributed subsequently to B. navacellensis n. sp. (from clade B) after COI sequencing; l-P, individuals 
sampled in the Folatiere cave and attributed subsequently to B. navacellensis n. sp. (from clade A) after COI sequencing; Q-V, individuals 
of an epigean population (Saint-Julien-de-la-Nef, Herault, France) of B. eurystoma. Scale bar: 1 mm. 


DISTRIBUTION. — The species is present in the subter¬ 
ranean aquifers of the Blandas and southern Larzac 
plateaus, France, Gard and Herault departments. Sample 
sites details are given in Table 1 and inferred distribution 
is shown in Figure 1. 

Description 

Shell 

Teleoconch conical and variably elongate, spire 
formed by c. 4 well convex whorls, gradually increas¬ 
ing, separated by deep sutures (Figs 4; 5). The shell 
width is c. % of the shell length and the last whorl is 
more than % of the shell length. Shell is transparent, 
with marked growth lines in the form of collabral 
grooves but no visible spiral micro-sculpture. Aperture 
holostome, oval; external margin of peristome slightly 


fluted. Umbilicus is widely open. Aperture could be 
rejected on the side, giving the shell a distinct globular 
shape. Operculum is thin, pliable, transparent and 
oval. Protoconch of \ l A whorls, almost smooth, but 
exhibits thin spiral micro-sculptures in the form of 
shallow grooves. Basic statistics of 11 shell parameters 
from 35 shells are given in Table 2. 

Anatomy 

The entire body is unpigmented. Eyes are invisible 
externally; however, a lens may be present but pig¬ 
mented retinal cells were not observed (Fig. 5C). Snout 
is about as long as wide, approximately parallel-sided 
or slightly tapered, with a weak but noticeable distal 
lobation. Tentacles are longer than the snout, more 
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Fig. 5. — Hypogean and epigean species of the Bythinella Moquin-Tandon, 1856 genus: A, holotype of B. navacellensis n. sp. (MNHN 
9315) from the Folatiere cave; B, live specimen of B. navacellensis n. sp. in situ from the type locality, with scattered fecal pellets visible 
on the ground and a drawn Buxus sempervirens L. leaf; C, histological section of the eye of a specimen from the Folatiere cave; D, shell 
of Bythinella eurystoma (Paladilhe, 1870); E, syntype of Bythinella padiraci Locard, 1903 (MNHN); F, holotype of Bythinella bouloti 
Girardi, Bichain & Wienin, 2002 (MNHN); G, holotype of Bythinella galerae Girardi, Bichain & Wienin, 2002 (MNHN). Abbreviations: 
L, lens; R, retina without pigmented layer. Scale bar: 1 mm. 


than 10 times as long as they are wide with the distal 
end slightly expanded. Under an optical microscope 
no ciliation is visible on the tentacles (Fig. 5B). 

On dissection, the anatomy of the male genitalia 
shows the typical organization found in Bythinella : 
penis bi-armed, the left arm (the penis proper) cy¬ 
lindrical and the right arm (usually referred to as the 
flagellum) short, blunt and thinner. The intestine 
makes a loop visible through the translucent shell. 
Content of stomach and intestine consists of clay 
and sand particles. 

Ecology 

The only known living population comes from 
the Folatiere cave. Living animals were found in 
shallow water, in streams and puddle in the active 
hydro-geological karstic network. In the Folatiere 


cave, about 67% of the individuals were found by 
sight trapped in puddle and 33% were collected 
on the submerged gallery walls, from the surface to 
about the depth of one meter. Only empty shells 
were collected by divers in deep water, i.e. more 
than one meter below the water surface. Since this 
species, which exhibits troglomorphic features, has 
only been sampled alive in hypogean habitats, we 
consider B. navacellensis n. sp. as a strictly stygobite 
species. 

DISCUSSION 

Introgression by hybridization 
According to Haase etal. (2007), genetic variability 
is usually low within the rissooidean populations and 
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specimens belonging to a single population cluster 
in the same clade after phylogenetic reconstruc¬ 
tion. Our phylogenetic results, based on the COI 
mitochondrial gene, highlighted that the hypogean 
population of the Folatiere cave was distributed across 
two well-supported clades separated by c. 4.7% of 
genetic divergence. A value higher than the level 
of inter-specific divergence (from 0.5% to 3.4%) 
commonly admitted for hydrobioids (Hershler et 
al. 1999, 2003; Liu et al. 2003; Hurt 2004). In 
contrast, phylogeny based on the ITS1 nuclear gene 
clustered the Folatiere cave sequenced specimens 
in the same clade. 

The topology conflict between nDNA and mtDNA 
trees could be explained by a scenario of introgres- 
sion by hybridization, i.e. a foreign genetic material 
transfer by interspecific backcross (Arnold 1997). 
Mitochondrial genome is considered as particu¬ 
larly susceptible to introgression events (Ballard & 
Whitlock 2004). This hypothesis implies that the 
misleading phylogenetic tree is given by the mito¬ 
chondrial gene(s), and that the “true” species tree 
is given by the nuclear genes. In our analyses, the 
hypothesis of introgressive hybridization involves 
only one evolutionary event corresponding to the 
mtDNA transfer from “ B . eurystoma' to “ B . nava- 
cellensis n. sp.” Following this scenario, the nuclear 
phylogenetic tree - the “right” evolutionary sig¬ 
nal - suggested that the Folatiere cave population 
constitutes a single and divergent evolutionary 
lineage. Shell characters, anatomic features (eye 
regression and unpigmented body) and ecological 
traits (cave habitat) corroborated this hypothesis 
as well. However, we admit that our scenario of 
introgression is based on an unbalanced molecular 
dataset. In order to test this hypothesis, new ITS1 
sequences should be produced to complement the 
current dataset. 

Distinction from other hypogean 
Bythinella species 

Bythinella navacellensis n. sp. is readily distinguishable 
on both morphological and molecular characters from 
the stygobite species B. padiraci (Fig. 5E) from the 
Padirac karstic network and morphologically from 
B. bouloti and B. galerae (Fig. 5F, G), all hypogean 
taxa belonging to the Atlantic drainage. 


Concerning the putative stygobite species B. eu- 
trepha , Bernasconi (2000) applied this name to 
epigean populations from southeastern France, but 
Falkner etal. (2002) restricted the name B. eutrepha 
to the hypogean population from the Lez basin. 
The latter considered that the material stored in 
the Bourguignat collection (Museum d’Histoire 
naturelle de Geneve) is composite with a specimen 
attributable to B. eurystoma (from St-Guilhem-le- 
Desert, Herault) (Fig. 5D) and a specimen cor¬ 
responding to the “true” B. eutrepha (from Lez 
alluviums). Consequently, Falkner et al. (2002) 
designated this specimen as the lectotype of B. eu¬ 
trepha. However, this type specimen was collected 
in the Lez rivers’ alluvium, making it impossible to 
assign it to any precise habitat. This species has never 
been sampled alive and its taxonomic status is still 
unclear. One of us (VP) has checked the Bythinella 
specimens in the Paladilhe collection (University of 
Montpellier) and collected recent shells from the 
alluviums of the Lez. Shells labeled “B. eutrepha ” 
in the Paladilhe collection seem composite with 
specimens from the Herault basin attributable 
to B. eurystoma and nine specimens from the Lez 
basin exhibiting distinct shell features attributable 
to the distinct morph B. eutrepha as it can still be 
sampled at present in the Lez alluviums. Overall, 
this material is morphologically very different from 
the new species described here. 

Distribution and IUCN category 
Among the 60 subterranean localities sampled, 
21 sites revealed empty shells attributable to the 
B. navacellensis n. sp. morphotype. The morph 
“B. navacellensis ’ is found in the karstic networks 
of the Herault basin, only in the karstic aquifers of 
the Causse de Blandas (Gard) and the Causse du 
Larzac meridional (Herault). It has not been found 
at altitudes below 300 m a.s.l. It seems to be absent 
from the Bueges river’s catchment and the upstream 
waters of the Herault catchment. No specimen was 
found in the hyporheic zone, even near the caves in 
which B. navacellensis vs. sp. occurs (e.g., in front of 
the Folatiere cave). Its distribution seems therefore 
to be limited by 1) the Atlantic drainage on the 
north-west side of its range; 2) the acid substrate of 
the Cevennes mountains in the north; 3) the Herault 
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river and the karstic network of the Bueges river in 
the East; and 4) lower altitude regions in the south, 
which are hydrogeologically isolated by a layer of 
marl. The infered distribution range (Fig. 1C) is 
a small delimited and cohesive area of c. 300 km 2 
located around the southern Larzac. 

Hydrological systems can be viewed as continental 
islands made of a single water body, isolated from 
adjacent ones. According to Prie (2008), within a 
single water body snails appear to be evenly distrib¬ 
uted in space and time in the whole hydrological 
system. Therefore, the distribution of subterranean 
snails (1) includes the whole water body they are 
found in and (2) can be delimited by adjacent 
isolated systems. The distribution map shown in 
Figure 1C is based on these assumptions. 

Population size or trend are unknown but pre¬ 
sumably stable in this relatively intact environment. 
However, short-range B. navacellensis n. sp. can be 
considered to occur in three locations only sensu 
IUCN (2001). According to the IUCN guide¬ 
lines, a location is defined as “a geographically or 
ecologically distinct area in which a single event 
(e.g., pollution) will soon affect all individuals of 
the taxon present”. Given the regional hydrologi¬ 
cal context, we reckon that three distinct locations 
harbour the whole species populations: 1) the 
Blandas plateau water body, which cohesiveness has 
been shown by fluorescein dye (Caumont 1992); 
and the 2) northward and 3) southward drainages 
of the Larzac plateau, which may constitute two 
separated water bodies. Bythinella navacellensis 
n. sp. can therefore be categorized as “Vulnerable, 
VU D2” (few locations) according to the IUCN 
criteria (IUCN 2001). 

Human activities are currently very limited in this 
area, with c. 1 inhabitant per km 2 , the main activity 
being extensive cattle farming. The aquifer is thus 
relatively free of pollution from human activities. 
Nevertheless, human pressure on southern France 
aquifers for drinking water is becoming a promi¬ 
nent issue and we cannot predict the consequences 
of such increasing pressure on either the aquifers or 
the stygobite fauna. This is a potential threat for 10 
species of stygobite snails found in the area (Prie un¬ 
published data), all of which being endemic to France 
and four of them endemic to the study area. 
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